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We calculate the Neel temperature Tjv for two-dimensional isotropic dipolar 
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expansion, employing the method of Callen. The theoretical predictions for for 
K 2 MnF 4 , Rb 2 MnF 4 , Rb 2 MnCl 4 and (CH 3 NH 3 )2MnCl4 are in good agreement with 
the measured values. 
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Recently |T| it has been shown that long-range order is possible in two-dimensional 
isotropic Heisenberg antiferromagnets due to the anisotropy of the dipole-dipole interaction. 
The occurence of a finite energy gap goes in hand with a nonvanishing order parameter and 
a finite Neel-temperature. In the present paper a quantitative improvement of the theory is 
achieved by means of Callen's extension of the Tyablikov decoupling scheme 0. 

In H we have used linear spin wave theory based on the Holstein-Primakoff transforma- 
tion to evaluate the magnon dispersion relation. The evaluation of the Neel temperature by 
means of the temperature independent dispersion relation leads to an overestimate of Tjy. 
In reality the magnon frequency softens with increasing temperature and thus the actual 
transition temperature is lower. This feature is accounted for by an extension of the Tyab- 
likov decoupling scheme due to Callen 0. In essence the dependence on the magnitude of 
the spin S is replaced by a, the temperature dependent order parameter. Even at T = the 
zero point fluctuations lead to a reduction of a as compared to S. The resulting transition 
temperature is lowered in comparison to the estimate of Ref. QXJ] such that a satisfactory 
agreement between theory and experiment is achieved. 

The Hamiltonian of a dipolar antiferromagnet reads 

l+V a/3 

with spins S/ at lattice sites x/. The first term in brackets is the exchange interaction Ju> 
and the second the usual dipole-dipole interaction. We consider a square lattice in the xy 
plane with lattice constant a and the spins orientated alternatingly along the z axis. The 
out-of-plane orientation is the classical ground state for the isotropic dipolar antiferromagnet 
with a nearest-neighbor exchange energy \J\ much larger than the dipole energy [|l|. 
Let us introduce the retarded double-time Green functions according to Callen || 

G\R k - R ,t) = ^Q(t)(lS£(t),e bS °S (0)]) = ((St\e bS °S^)), (2a) 
G 2 (R,-R ,t) = ^Q(t)(iS k (t),e bS °S (0)]) = ((S k \e bS °S )), (2b) 

which obey the following Fourier transformed equation of motion |||4j 
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^m\^So)) u = e 6 + ((lS+,H]\e bS SS }}^ , (3a) 
u((S^S^)) u = m,H]\e bS o S ^ , (3b) 

with the equal time commutator 

Q b =([S\e bSZ S-]). 

Here R is a lattice point of sublattice one and is a general lattice vector. The higher 
order Green functions generated by the commutator in Eqs. (3a) and (3b) are approximated 
by the Tyablikov decoupling scheme 

mS+leWSo)) - (Sf)(<S+|e 6S °S ->>- (4) 

As a consequence of translational symmetry the mean spin value is independent of the lattice 
site for each sublattice: (SQ = —{SQ = a where li E L x (l 2 £ L 2 ) refers to sublattice one 
(two). Using this approximation we obtain from Eqs. (3a) and Eq. (3b) a set of four 
equations for the Green functions G q , Gq +qo , G q and Gq +qo . The evaluation of the order 
parameter and the magnon dispersion relation requires only two of them: 

GqH = + C? q+q „H) = b (^2_ e 2 )(cL , 2 _ e 2) > ( 5 ) 



with 



" - 2^q^q+q« ~~ ^q+qo) + 2^+1° ^q ~~ ^q) 

+K + cj +qo )(A q + A q+qo ) + (cj - ^ +qo )(s q+qo - s q ), 

= ~~ ^ + ^q+qo ~ B l+10 + SCqCq+qo)) 

C (1) = i(A q + A q+qo ). 

Here q denotes the wave- vector of the chemical Brillouin zone and q = ^(1,1,0). In 
the magnetic Brillouin zone which is half the chemical there exist two distinct spin-wave 
branches with frequencies (i — 1,2) which read 

4 = l^i ±^2), (6) 



^1 = B\ + ^4q+ qo _ ^q+q + ^Cq ^q+qo > 

^2 = (^q — ^q — ^q+qo + "^q+qo) 2 + [C q+qo (y4 q+q() — -B q+qo ) — Cq (A q — Bq)} 
X [Cq (^4q+q + -^q+qo) — C q+qo (Aq + Bq)]. 

In Eq. (6) the coefficients 

Aq = cr(2J qo — Jq — Jq+q ) + a(2A^ Q — A^ — A q + qo ), (7a) 
Bq = a(Jq + q - Jq) + a(A^ - A™), (7b) 
Cq = iaA x «, (7c) 

have been introduced. This result, Eq. (6), coincides with the magnon frequencies derived 
by the Holstein-Primakoff transformation when a is replaced by S and the external magnetic 
field is set to zero But there is a difference even at absolute zero temperature, because 
the ground state of the antiferromagnet is not the Neel state of fully aligned spins, i.e. 
a(T = 0) < S as will be seen later. Note that the magnon frequency scales with the order 
parameter a, i.e. the whole spectrum softens with increasing temperatures and vanishes at 
the phase transition. 

Now we turn to the evaluation of a. For arbitrary spin S the spin expectation value is 
given by the well known relation 

(S Z )=S(S+1)-((S Z ) 2 )-(SS + ). (8) 

For 5=1/2 the order parameter can be calculated directly via the Green functions of Eq. 
(2a) and (2b) with 6 = and Oo = 2a. For higher spin quantum numbers the above formula 
is not so helpful. Then a convenient starting point is the following gerneralized thermal 
average 

m = (e bS »S S+), (9) 
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by means of which a self-consistent system of equations can be derived by the method of 
Callen. The above thermal average can be represented by the spectral theorem @,|]: 

(e^U)=lnn-- J du, n(u)- £{Gj( W + iS) - G\(u - i5)} 



= e 6 (n-l/2), (10) 

where Eq. (5) has been used in the last step and N denotes the total number of spins. Here 
we have introduced the Bose occupation number 



n 



and 



» = jfX M^j + 7M=4) " (ei,+ ^-«!) " (<2) - (U) 

The right hand side of Eq. (11) depends on o only via the occupation numbers n(ej). For 
spin 5 = 1/2 and vanishing parameter 6 the thermal average [Eq. (9)] represents the number 
of spin wave excitations, which reduce the staggered magnetization from the totally ordered 
Neel state, not only for finite but also for zero temperature. 

For arbitrary spin one has to express a in terms of n which can be achieved by the 
method of Callen with the result : 

1 (n + l/^ 1 + (n - 1/2)^ 
° ~ {b + 2 j (n + 1/2)^+1 _ ( n _ TJ^S+i n ' ^ 

Eqs. (11) and (12) constitute a self-consistent system of equations for n and the spin 
expectation value a. 

Let us now discuss the dispersion relation for T = 0. In this limit Eq. (11) reduces to 

- it m 1 - All> ± C '"" 62 
no = n(T = )=-^ 2eif2(ei+£j) , (13) 

where the right-hand side is independent of a. Let us denote the spin expectation value at 
zero temperature by cr = cr(T = 0) which is found by inserting Eq. (13) into (12). Knowing 
ctq one obtains for the staggered magnetization N(0) at T = 0: 



N(0) = gfi B Na . (14) 

One can convince oneself that Eq. (14) for large S coincides with the expression derived by 
the Holstein-Primakoff transformation [TLJ (see also U). This must be so because the latter 
is an expansion in 1/S. In addition we derive an energy gap (q = 0) from Eq. (6) 



Eg ee ei (T = 0) = e 2 (T = 0) = 2a oy /A« - A q ^(J qo - J ) - « - A") , (15) 

with 

APP — A^x _ AW APP — A xx — AW 

^0 — ^0 — ^0 5 ^qo — ^qo — ^qo- 

This is of the same form as the result from the Holstein-Primakoff transformation Eq |I[ 
except for the prefactor which is smaller by the ratio ar /S. 

Now we turn to the evaluation of the transition temperature T/v, i.e. consider the limit 
a — > 0. Since the spin-wave energy (Eq. 6) is proportional to o the Bose occupation numbers 
can be replaced by their classical limit: 

/ \ k B T 
n[€i) -> . 

If this is inserted into Eq. (11) together with o —>■ one obtains 

k B T N -AV 
n =^J-Y.7TJ^ ■ ( 16 ) 



q y-i'-tj 

To keep track of the a dependence we have introduced the cx-independent quantities A^ 1 ' = 
A^'/a 3 and q = e^/a. According to Eq. (16) n increases indefinitly with a — > and thus 
the second relation between a and n, Eq. (12), becomes 

,= *£±ili. (17) 
3 n 

Combining Eqs. (16) and (17) we obtain an explicit expression for the Neel temperature: 

3k B 

with 



For purely isotropic antiferromagnets the coefficient F diverges, excluding long-range or- 
der at finite temperature in two dimensions in accord with the Hohenberg-Mermin- Wagner 
theorem ||[7|. 

In the presence of the dipolar interaction there is an energy gap and TV becomes finite. 
If the dipolar interaction is weak in comparison with the exchange energy {g^B) 2 /o? "C \ J\ 
and if the argument of the summation is approximated by its small q limit one obtains 

T n ~ 7^ , (19) 

which coincides with the analogous formula derived by the Holstein-Primakoff transforma- 
tion (Eq. (18) in Ref. JIJ). In the general case the above sum [Eq. (18)] is evaluated with 
the full dispersion relation Eq. (6) and by computing 100 x 100 points in the 2D Brillouin 
zone and determining the other points by linear extrapolation. The dipole sums have been 
calculated via Ewald summation ||. 

Now we apply our theory to real quasi two-dimensional antiferromagnets. Prominent 
examples of almost two-dimensional antiferromagnets are the tetragonal antiferromagnetic 
halides K 2 MnF 4 , Rb 2 MnF 4 , Rb 2 MnCl 4 and (CH 3 NH 3 ) 2 MnCl4. In these quadratic layer 
structures the out-of-plane exchange interaction is neglegible in comparison to the in-plane 
exchange interaction (about 10~ 4 ) |],|nj whereas the dipole energy is larger by an order of 



magnitude. This two-dimensional character has been shown experimentally by the absence 
of any dispersion along the z-direction [ 1(J . For these halides the measured exchange energy 



| J\, the lattice constant a, the energy gap Eq XP , the spin-flop field H^ p and the transition 
temperature T^ xp [|10| -|T2"|| are listed in table I. The spin-flop field H s f is the critical magnetic 



field at which the antiferromagnetic Neel ground state changes to the spin-flop ground state. 



It can be calculated by adding to the Hamiltonian, Eq. (1), the Zeeman energy [|T^]. From 
the full dispersion-relation the spin-flop field is defined by that field for which the magnon 
energy vanishes 



H « = 7^r E °> (20) 

in close analogy to the formula obtained by the Holstein-Primakoff method (Eq. (11) in 
Ref. The table also contains the theoretical energy gaps Eq and Eq calculated via 

Eq. (15) and Eq. (9) of Ref. 0, the resulting spin-flop fields H£ and H sf via Eq. (20) 
and Eq. (11) of Ref. [0J and the theoretical transition temperature Tjv [Eq. (18)]. All 
these substances have spin S = 5/2, which yields from Eq. (13) and (12) o"o = 2.30 for pure 
isotropic antiferromagnets. This value is increased only neglegibly by the dipolar interaction 
as can be seen from table I. 

We find a good agreement with the measured the Neel temperature although our theory 
accounts only for the dipolar interaction and no other anisotropy. Corrections due to dipolar 
interactions between different planes are neglegible because of the large lattice constant in 
z-direction; e.g. the energy gap for K 2 Mn0 4 [Eq. (6)] is altered by 

E° (3D) = E°q(2D)[1 + O(10-% (21) 

which justifies the application of the two-dimensional model. Note that the nearest-neighbor 
exchange energy J is the only parameter entering in our theory. Experimentally this param- 
eter has been derived by fitting the measured spin wave spectrum with a dispersion relation 
which is different of ours [Eq. (6)]. 

For the halides listed in table I, the energy gap obtained from Eq. (15) and the transition 
temperature are lower than the experimental values. The following reasons may be respon- 
sible for that: (i) In the Holstein-Primakoff approximation the softening of the magnons is 
neglected entirely. This leads to an overestimate of T^. In the Callen method the magnons 
soften in the entire Brillouin zone, thus particularly near the phase transition the softening 
is overestimated and leads to a which is somewhat too low. (ii) A small readjustment of 
J could be necessary if our dispersion relation, Eq. (6), is used to fit the data, (iii) A small 
additional anisotropy from the crystal field might be present as suggested by [|HH . 



8 



ACKNOWLEDGMENTS 

This work has been supported by the German Federal Ministry for Research and Tech- 
nology (BMFT) under the contract number 03-SC2TUM. 



9 



TABLES 

TABLE I. Exchange energy |J|, lattice constant a, energy gap £"0, spin-flop field H s {, Neel 
temperature Tjv and zero temperature order parameter ctq. 





\J\ 


a 




E o 




sf 


H J 




rpexp 

N 


rpth 

N 


O"0 




[K] 


[A] 


[K] 


[K] 


[K] 


[T] 


[T] 


[T] 


[K] 


[K] 




K 2 MnF 4 


8.5 a 


4.17 a 


7.4 b 


7.1 


7.6 


5.4 d 


5.3 


5.6 


42 a 


41 


2.33 


Rb 2 MnF 4 


7.4 C 


4.20 g 


7.3 b 


6.5 


7.0 




4.8 


5.1 


38^ 


36 


2.33 


Rb 2 MnCl A 


11. 2 f 


5.05 e 


7.5 f 


6.1 


6.6 




4.5 


4.9 


56 f 


48 


2.32 


{CH 3 NH 3 ) 2 MnCl 4 


9.0 f 


5.13 c 




5.3 


5.7 




3.9 


4.3 


45 f 


39 


2.32 



a Reference |1 
b Reference |lj| 



c Reference || 
rf Reference [11] 
e Reference 
f Reference (IT 



9 Reference [17] 
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